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(57) Abstract: The invention concerns a method for processing a signal using an approximate MAP (maximum a posteriori) al- 
gorithm for determining a likelihood ratio A^* of a set of states X of a lattice at a time k, with each of said states being associated 
at least one intermediate variable belonging to a group comprising a so-called forward variable and a so-called backward variable, 
propagated by said MAP algorithm and recursively calculated respectively in a direct orientation and in an indirect orientation at 
said time k relative to said lattice. The invention is characterized in that said process comprises a step which consists in reducing the 
numb^ of selected states by said MAP algorithm so as to calculate said likelihood ratio, and, for at least some unselected states, in 
assigning to said forward variable and/or said backward variable at least one specific value, to calculate an approxiD:iate likelihood 
ratio. 
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(57) Abr^g^ : Proc^6 de traitement d'un signal mettant en oeuvre un algorithme de type MAP approch6 et applications corres- 
pondantes. Uinvention conceme un proc6de de traitement d*un signal mettant en oeuvre un algorithme de type MAP ("Maximum a 
posteriori") permettant de determiner un rapport de vraisemblance A^^ d'un ensemble d'^tats X d'un tieillis k un instant k, h. chacun 
desdits 6tats 6tant associ^es au moins une variable interm^diaire appartenant au groupe comprenant une variable appel^e "fonvard" 
et une variable appelde "backward". propag6es par ledit algorithme MAP et calcul^es r^cursivement respectivement dans un sens 
direct et dans un sens indirect audit instant k par rapport audit treillis- Selon Tinvention, un tel proc6d6 comprend une 6tape de 
reduction du nombre d'6tats selectionn^s par ledit algorithme de type MAP en vue d*un calcul dudit rapport de vraisemblance, et, 
pour au moins certains 6tats non-s61ectionn6s, on affecte k ladite variable "forward" el/ou "backward" coirespondante au moins une 
valeur d6termin&, de fagon k calculer un rapport de vraisemblance approch^. 
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Procdd^ de traitement d'un signal mettant en oeuvre un algorithme de 
type MAP approch^ et applications correspondantes. 

Le domaine de rinvention est celui du traitement de signaux transmis dans 

un canal perturb6, ou susceptible de Tdtre. 
5 Plus pr6cis6ment, rinvention conceme un proc6d6 de traitement du signal, 

mettant en ceuvre un algorithme de type MAP (pour "Maximum A Posteriori"). 

Ualgorithme MAP est notanmient pr6sent6 dans Tarticle de L, E. Baum et T. 

Pertie, intitule "Statistical Inference for Probabilistic Functions of Finite State 

Markov Chains" (en frangais, "Interference statistique pour les fonctions 
10 probabilistes de chalnes de Markov k 6tats finis"), Ann. Math. Stat. 37 : 1554- 

1563, 1966. 

L'algorithme MAP, aussi appel6 algorithme du Forward-Backward (FB), 
suscite actuellement beaucoup d'int6r8t. En effet, cet algorithme permet 
notamment de decoder un code convolutif en associant une information de 

15 fiabilite aux bits decodes, et de nombreuses applications, telles que celles mettant 
en ceuvre un "turbo-code" par exemple, peuvent en titer profit 

On rappelle rapidement ci-apr6s le principe de Talgorithme MAP. Par 
algorithme MAP, on entend ici, et dans toute la suite du document, aussi bien les 
algorithme de type MAP, que les algoridmies qui en sont d6riv6s, et notamment 

20 les algorithmes de type Log-MAP et Max-Log-MAP. 

De manifere sch6matique, I'algorithme MAP pent Stre consid6r6 comme la 
fusion des informations Qes m6triques de noeuds) g6n6r6es par deux d6codeurs de 
Viterbi modifies travaillant en sens inverse sur les donn^es. On rappelle que 
Talgorithme de Viterbi, couramment utilise dans le domaine du d6codage de 

25 signaux regus, consiste h, rechercher dans un treillis repr^sentatif du codage, le 
chemin qui correspond ^ la sequence d*6tats la plus probable, c'est-£t-dire ceUe qui 
est k la distance minimale de la sequence regue. Cet algorithme est d6crit dans 
Tarticle de G. D. Forney, "The Viterbi Algorithm", Proceedings of the TRRR, vol. 
61, No. 3, mars 1973. 



wo 03/075468 PCT/FR03/00679 



L'algorithme MAP pennet de daenniner le rapport de vraisemblance de la 
prpbabiUt^ a posteriori des 6tats k r instant k, 0 ^ A ^ , d'une chaine de Markov 
cach6e k 6ms finis, k partir d'observations Yf =^„,...,Yj^y On d6finit 

g6h6ralement una chaine de Markov comme une collection d*6tats al6atoires 
5 cons^cutifs enchatoes les uns aux autres, dans laquelle le passage d'un 6tat k un 

autre est associ6 k une probability, appel6e probability de transition. 

On pent ainsi considerer une chaine de Markov {Xf^}k valeurs dans un 

espace fini, £ = ^,-JV^}, d^finie par une matrice de transition entre etats 

n = (;r^y), de dimension NxN et une loi initiale v = (v^X^i d6finie par 

10 Vf = Pvob[X^ = f], Vi e E. Selon la terminologie classiquement employee, une 

transition entre deux 6tats X^.! et Xj, est d^finie par un couple d'6tats, et est notfe 
= (A'^^i, JSr^) . La matrice 11 est la matrice de transition entre 6tats assocife k 

{Xf^} dont les Sl^ments Ttj j sont d6finis par : 

n^ j = ProblX^ = j/X^_^ = i], pour tout ij e E. (1) 

15 Dans toute la suite de ce document, on considfere le cas des modMes de 

Markov cach6s, pour lesquels on dispose d' observations bruit6es ^j,} k valeur 

dans SI"*, recueillies k travers un canal sans memoire. En d'autres termes, 
conditionnellement aux transitions entres 6tats 4f^, les observations {7^} sont 
mutuellement independantes, et chaque observation Yj^ ne depend que de la 
20 transition au mSme instant k. Par consequent, le module de signal observe pent 
s'^crire sous la forme suivante : 

Y,^Hi^,) + B, (2) 

oil H(.) repr6sente la fonction de transfert du canal sans m6moire d^fini dans ^ 
et k valeur dans 91'', = (JST^t^j, JT^) est un couple d'etats k r instant k k valeur 
25 dans ExE , 7^ repr6sente une observation k Tinstant k k valeur dans SH'' et ot Bj, 

est un bruit additif blanc. 

Dans le contexte des conmoiunications num^riques, le module d^fini par 
r^quation (2) peut s'appliquer en ^galisation : le canal ifQrepr6sente alors le 
canal de transmission en bande de base et Xj^ repr&ente une suite de symboles de 

30 r alphabet de modulation dont la longueur est ^gale k la memoire du canal. 
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Lors du d^codage d'un code convolutif, H{.) repr6sente alors les 
polyndmes g6n6rateurs du code, et Xf^ repr&ente le contenu de la mSmoire du 

codeur de canal. 

On d6crit d6sormais plus en details la notion de rapport de vraisemblance 

5 calculi dans le cadre de Talgorithme MAP, 

Le rapport de vraisemblance d'un ensemble -Yappartenant k Tensemble 

des parties de E est d6fini conune le rapport de la probability a posteriori de 
Tensemble X sur la probability a posteriori du compl6mentaire de X dans E. 

Le rapport de vraisemblance pent encore 6tre expiim^ sous la forme de 
10 r^quation suivante : 

oil PAP/ = P[Xj, = i/Y^ = yo»— ^y^],0<k<K . 

L'algorithme MAP est classiquement associ^ k une representation en 
treillis des diff6rents 6tats k consid6rer. Ainsi, si Ton considfere un codeur 

15 convolutif, Talgorithme MAP utilise, pour le d6codage, la representation en treillis 
de revolution temporelle du codeur. Le treillis est generalement repr6sent6 par un 
graphe, dont Taxe des abscisses figure le temps et Faxe des ordonnees repr6sente 
les difiKrents etats possibles. Chaque noeud du treillis est done associe h un etat 
d'une part, et 2t un instant k d'autre part Les transitions possibles entre les noeuds 

20 de deux instants successifs k et k+1 sont repr^sentes par des branches. 

L'algorithme MAP fait intervenir deux variables interm6diaires, 
respectivement appel^es variable "forward" et variable "backward". Ces variables 
propagees par Talgorithme MAP sont calcul6es recursivement, respectivement 
dans le sens direct ("forward") et indiiect ("backward") par rapport au treillis, et 

25 permettent d'aboutir au calcul du rapport de vraisemblance k chaque instant k . 

Ces variables "forward" et '"backward", ddjOnies dans la suite du document, 
utilisent la probabilited' observer y^^ conditionnellement aux 6tats Xk=j et Xjt-i=f, 

exprim^e par I'^quation suivante : 

b,j(X,) = VTob[Y, = y, /X,,, = hX, = J] (5) 
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Cette probability est classiquement appel6e mStrique. Elle peut 6ventuellement 
ddpendie de k» et est fonctioii de la density de probability du bruit B^^ . 

Ualgorithme MAP procfede en trois 6tapes, h savoir une 6tape de recursion 

avant (ou "forward"), une etape de recursion arrifere (ou "backward") et enfin, une 

5 6tape de fusion des donnees. 

La recursion "avant" consiste h calculer les probabilit6s af de passer par 

les noeuds repr&entatifs des ^tats Xi=i aux instants k, k partir des informations 
obtenues en parcourant le treillis dans le sens direct La variable "forward" a* est 
ainsi d^finie conome la probability de ryvenement ( -ST, = /,7o = J'o'— = -V* ) • 

10 af = Pr ofe[J5r, = /.Fo = yo,-^>Y, = y,] (6) 

On peut d^montrer que les a*, 0^k<K ^sc calculent r6cursivement de 
la f agon suivante : 

avec la valeur initiale af = Vp Vi e E, 
15 De m6me, la recursion "arrifere" consiste i calculer les probabilit^s jS* de 

passer par les noeuds repr^sentatifs des ^tats Xi=i aux instants k, k partir des 
informations obtenues en parcourant le treillis dans le sens indirect. La variable 
"backward" j3* est ainsi dyfinie comme la probability de ryvynement 

20 j8f = PxoblX, = UY,,, = = VkI (7) 

On peut dymontrer que les j3 * , 0<^k<K calculent rycursivement de 

la f agon suivante : 

avec la valeur initiale jSf = 1, Vi e J?. 
25 La troisi^me ytape mise en oeuvre par Talgorithme MAP consiste k 

fusionner les rysultats des deux rycursions pour calculer la probability "a 
posteriori" des ytats, en calculant le rapport de vraisemblance des symboles 
dytectys. 
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Le rapport de vraisemblance d'un ensemble d'6tats d6fiiii par I'equation (3) 
peut alors 6tre exprim6 en fonction des variables interm^diaires "forward" et 
"backward" sous la forme : 

5 Dans le contexte de la detection de symboles par exemple, Tensemble X 

regroupe les Stats correspondant h remission d'un m8me symbole de la 
modulation. Dans le contexte de d^codage de canal, Tensemble X regroupe les 
6tats correspondant au mgme mot a I'entr^e du codeur de canal. 

L'algorithme MAP, bien que foumissant des r6sultats tr&s performants, 
10 pr6sente rinconv6nient d'Stre trfes complexe k mettce en oeuvre. 

Notamment, la realisation mat^rielle de l'algorithme MAP est trfes delicate, 
et les r&epteurs impl6mentant cet algorithme sont par consequent coflteux. 

Les problemes d'impldmentation de Talgorithnie MAP classique sont en 
partie dus k la forte dispersion des probabilites * qui n^cessite done une grande 
15 precision pour le calcul des variables "forward" a* et "backward" jS * . En outre, 

le calcul du rapport de vraisemblance tel que d6fini par I'equation (3) ci-dessus 
impose d'effectuer des divisions qui sont cofiteuses et peuvent poser des 
probl&mes de stabilit6. 

Pour resoudre ces differents inconvenients de l'algorithme MAP, on a 
20 propose des algorithmes de types Log-MAP ou Max-Log-MAP. 

Neanmoins, les algorithmes de type Log-MAP ou Max-Log-MAP, comme 
les algorithmes de type MAP, restent txhs complexes a mettre en ceuvre, et 
necessitent d'effectuer un grand nombre d'operations lors du calcul du rapport de 
vraisemblance. 

25 L'invention a notamment pour objectif de pallier ces inconvenients de I'art 

ant^eur. 

Plus predsement, un objectif de Tinvention est de foumir une technique de 
traitement de signal mettant en oeuvre un algorithme MAP k complexite reduite 
par rapport aux techniques de I'art anterieur. Notamment, l'invention a pour 
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objectif de permettre une reduction du nombre de calculs mis en cBUvre, taut par 
les algorithmes de type MAP, que par les algorithmes de type Log-MAP ou Max- 
Log-MAP, en r^alisant un meilleur rapport performance / complexit6 que dans 
Tart ant6rieur. 

5 Un autre objectif de I'invention est de mettre en oeuvre une telle technique, 

permettant la conception de r^cepteurs moins complexes, et done moins coiiteux 
que dans Tart ant^rieur. 

L'invention a encore pour objectif de foumir une telle technique qui 
realise un compromis satisfaisant entre les exigences de faible complexity d*une 

10 part et de fiabilit6 des r^sultats d'autie part. Ea d*autres tennes, Tinvention cherche 
k foumir une m6thode permettant de r^duire syst6matiquement, avec une bonne 
quality d'approximation, la complexity des algorithmes de types MAP, Log-MAP 
et Max-Log-MAP. 

Ainsi, alors que les algorithmes de type MAP classiques n^cessitent 

15 g6neralement la mise en ceuvre d'un nombre d*op6rations en A^^ oil L repr^sente 
la longueur du canal de transmission, et oii A repr^sente le nombre de caractdres 
de Talphabet de modulation, Tinvention vise k foumir une technique mettant en 
ceuvre un algorithme de type MAP mettant en ceuvre un nombre d'op^rations en 
A^, oil M<L-1, tout en conservant une bonne quality d'approximation du rapport 

20 de vraisemblance. 

L'invention a ygalement pour objectif de mettre en ceuvre une telle 
technique qui permette une prise de dycisions souples. 

Ces objectifs, ainsi que d*autres qui apparaitront par la suite, sont atteints h 
Taide d'un procydy de traitement d*un signal mettant en oeuvre un algorithme de 

25 type MAP ("Maximum a posteriori") permettant de dyterminer un rapport de 
vraisemblance A^^ d'un ensemble d'ytats X d'un treillis k un instant k, k chacun 
desdits ytats ytant associyes au moins deux variables intermydiaires, 
respectivement appeiyes variable "forward" et variable "backward", propagyes par 
ledit algorithme MAP et calcuiyes rycursivement respectivement dans un sens 

30 dkect et dans un sens indirect audit instant k par rapport audit treiUis. 
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Selon rinvention, un tel proc6d6 comprend une 6tape de reduction du 
nombre d'6tats s61ectioim6s par ledit algoritfame de type MAP lors d'un calcul 
dudit rapport de vraisemblance, et, pour au moins certains dtats noii-s61ectioim6s, 
on affecte k ladite variable "forward" et/ou "backward" correspondante au moins 

5 une valeur d6termin6e, de fagon h calculer un rapport de vraisemblance approch6. 

Ainsi, rinvention repose sur une approche tout k fait nouvelle et inventive 
de mise en oeuvre d'un algorithme MAP, pemiettant de simplifier les diffdrentes 
Stapes de traitement des signaux tout en maintenant une quality de r6sultats 
satisfaisante. En effet, rinvention propose notamment de reduire le nombre d*6tats 

10 s61ectionn6s par Talgorithme MAP pour le calcul du rapport de vraisemblance, de 
fagon k simplifier les diff6rentes Stapes mises en oeuvre lors d'un tel calcul, et 
d'attribuer aux variables intermediaires manquantes des valeurs d6termin£es 
pemiettant de calculer un rapport de vraisemblance approch6. L'invention repose 
done sur un algorithme de type MAP approch6, pr6sentant un coiit de calcul r6duit 

15 par rapport aux algorithmes MAP classiques. Les recepteurs de signaux mettant 
en oeuvre un tel algorithme sont done moins complexes k construire que les 
recepteurs de Tart ant6rieur. 



Avantageusement, pour un instant k donn6, ladite au moins une valeur 
determin6e a(k) affect^e k ladite variable "forward" est telle que 




ladite variable "backward" est telle que 0^b(k)^Min(P^), ot et M^^ 



representent un ensemble desdits 6tats s61ectionn6s respectivement dans ledit sens 
direct et dans ledit sens indirect audit instant k, et oil al" et repr6sentent 
respectivement lesdites variables "forward" et "backward" audit instant k. 
25 Pr6f6rentiellement, pour un instant k donn6, ladite valeur detennin6e a(k) 

et/ou b(k) est unique et est affect6e k au moins une variable "forward" et/ou 
"backward" ft^ 

Selon une variante avantageuse de rinvention, on affecte une valeur 
constante auxdites variables "forward", respectivement "backward", de fagon que 
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ledit algorithme de type MAP soit un algorithme unidirectionnel de type direct, 
respectivement indirect 

Selon une caract6ristique avantageuse de rinvention, ladite 6tape de 
reduction du nombre d'^tats met en oeuvre un algorithme de recherche en treillis 
5 de type "breadth-first". 

Selon une premiere variante pr6f6rentielle, ledit algorithme de type 
"breadth-first" est un algorithme de M, 

Selon une deuxifeme variante pr6f6rentielle, ledit algorithme de type 
"breadth-first" est un algorithme de type T, utilisant au moins un seuil. 
10 Pr6f6rentiellement, ledit au moins un seuil est variable en fonction dudit 

instant k. 

Avantageusement, on affecte audit seuil variable une valeur pr6d6tennin6e 
pour chaque instant k. 

De mani^re pr6f6rentielle, pour chaque instant k, la valeur dudit seuil 
15 variable est determin6e par mise en oeuvre d'un algorithme adaptatif . 

De manifeie avantageuse, ledit algorithme adaptatif est du type algorithme 
de gradient. 

Selon une caract6ristique avantageuse de Tinvention, ledit treillis 
comprenant une pluralit6 de noeuds associ6s chacun k un desdits 6tats et ^ un 
20 instant k donn6, la valeur dudit seuil variable T k un instant 0^+1) est d6tennin6e 
par r^quation suivante : 

T(k + 1) = Tik) - fi{Mik) - M,) 
ot T(k) repr^sente la valeur dudit seuil variable audit instant k, 
Mc est le nombre de noeuds dudit treillis propag6s visd, 
25 M(k) est le nombre de noeuds dudit treillis propag6s audit instant k, 

et \i est une constante positive representative d'un gain d'apprentissage. 
Preferentiellement, ledit algorithme adaptatif est du type algorithme de 

gradient h, pas variable. 

Avantageusement, ledit gain d'appientissage \i est fonction dudit instant k. 
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Pr6f6rentieUement, ledit algorithme de type "breadth-first" 6tant un 

algorithme de type M, lesdites valeurs d6temii]i6es a(k) et/ou bCk) affect^es 

respectivement auxdites variables "forward" et/ou "backward", k un instant k 

domi6 sont denudes par les Equations suivantes : 
5 aik) = Min(.a^)-c. 

iGMt 

ot Cf et C5 sont deux constantes positives. 

De manidre pr6f6rentieUe, ledit algorithme de type "breadth-first" 6tant un 
algorithme de type T, lesdites valeurs d6termin6es a(k) et/ou b(k) affectfes 
10 respectivement auxdites variables "forward" et/ou "backward", k un instant k 
donn6, sont doim6es par les ^nations suivantes : 
a{k)^T^ik)-Cf 
b(k) = T\k)-c, 

oil Cf et Ci, sont deux constantes positives, et oil T^(k) et T^(k) d6signent la valeur 
15 dudit seuil variable audit instant k respectivement dans ledit sens direct et dans 
ledit sens indirect 

Avantageusement, ledit algorithme de type MAP appardent au groupe 
comprenant : 

les algorithmes de type MAP ; 
20 - les algorithmes de type Log-MAP ; 

les algorithmes de type Max-Log-MAP. 

On notera que les variables "forward" et "backward" peuvent avoir des 
definitions diff6rentes selon le type d'algorithme utilis6 (MAP, Log-MAP ou 
Max-Log-MAP), ainsi que d6taiI16 dans la suite du document 
25 Pr6f6rentiellement, ledit algorithme de type MAP 6tant un algorithme 

unidirectionnel, ledit proc6d6 met en ceuvre une 6tape de comparaison de 
decisions prises par ledit algorithme unidirectionnel avec des decisions 
correspondantes prises par un algorithme de type Viterbi, appel6es d&isions de 
Viterbi. 
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Par algorithme de type Viterbi, on entend ici, et dans toute la suite du 
document, un algorithme de Viterbi k nombre d'6tats r6duit, auquel on associe une 
operation de "trace-back" (c'est-a-dire de d^codage de Talgorithme). Ainsi, les 
decisions de Viterbi sont des decisions foumies par un algorithme k nombre 
5 d'6tats r6duit, auquel est associe un "trace-back", 

Selon une premiere caract^ristique avantageuse, en cas de comparaison 
negative pour au moins une desdites decisions prises par ledit algorithme 
unidirectionnel, ledit proc6d6 met en oeuvre une 6tape de substitution de ladite 
decision de Viterbi correspondante k ladite decision prise par ledit algorithme 
10 unidirectionnel, appel^ d^ision substitu^. 

Pr6f6rentiellement, on affecte k la valeur absolue dudit rapport de 
vraisemblance associ6 k ladite decision substituee une valeur d6termin^ V. 

Avantageusement, ladite valeur d6termin6e V est 6gale k la valeur absolue 
du rapport de vraisemblance moyen de la sequence, 
15 V correspond done k la moyenne des fiabilit6s, soit environ k la fiabilit6 

moyenne de la sequence. 

Selon une deuxifeme caract6ristique avantageuse, en cas de comparaison 
negative pour au moins une desdites decisions prises par ledit algorithme 
unidirectionnel, ledit proc6d6 met en oeuvre une 6tape de pond^ration dudit 
20 rapport de vraisemblance associ6 k ladite decision consid6r6e, tenant conq>te de 
ladite decision de Viterbi. 

Pr^fSrentiellement, Y 6tant un ensemble d'6tats associ6s k une decision 
foumie par ledit algorithme de type Viterbi k un instant i, et repr6sentant le 
rapport de vraisemblance associ6 kY k Tinstant i, calcul6 par ledit algorithme 
25 unidirectionnel, au cours de ladite 6tape de ponderation, on remplace Ai^ par A| 
d^fini par A^ = A^ + Z3f x V, oil V est une valeur d6termin6e. 

^invention s'applique avantageusement k Tun au moins des domaines 
appartenant au groupe comprenant : 
la detection de symboles ; 
30 - le codage/ddcodage de signaux ; 
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le turbo-d^codage ; 
la turbo-d6tection ; 

le codage de source par quantification (notamment quantification en 
treiUis). 

5 L'invention concerne 6galement un r6cepteur de signaux de 

communication comprenant des moyens de mise en ceuvre d*un algoritiime de 
type MAP ("Maximum a posteriori") permettant de determiner un rapport de 
vraisemblance A|,^ d'un ensemble d*6tats X d'un treillis k un instant k, 
k chacun desdits Stats 6tant associ6e au moins une variable interm6diaire» 
10 appartenant au groupe comprenant une variable appelSe "forward" et une variable 
appelSe "backward", propagSes par ledit algorithme MAP et calcul6es 
r6cursivement respectivement dans un sens direct et dans un sens indirect audit 
instant k par rapport audit treillis. 

Selon rinvention, un tel r^cepteur comprend des moyens de reduction du 
15 nombre d*6tats s61ectionn6s par ledit algorithme de type MAP en vue d'un calcul 
dudit rapport de vraisemblance, et, pour au moins certains 6tats non-s61ectionn6s, 
on affecte k ladite variable "forward" et/ou "backward" correspondante au moins 
une valeur dStenninSe, de fa^on k calculer un rapport de vraisemblance approch6. 
D'autres caractSristiques et avantages de rinvention apparaitront plus 
20 clairement k la lecture de la description suivante d'un mode de realisation 
pr6f6rentiel, donnS k titre de simple exemple illustratif et non limitatif, et des 
dessins annexes, parmi lesquels : 

la figure 1 prSsente un treillis reprSsentatif d'une plurality d'Stats associS k 
un algorithme MAP de rinvention ; 
25 - la figure 2 prSsente un organigramme illustratif du principe de "trace-back" 
comparatif ou comparatif pondSrS, mis en ceuvre selon rinvention en 
combinaison avec les algorithmes MAP unidirectionnels ; 
la figure 3 pr6sente un synoptique des difi^6rentes Stapes de calcul mises en 
oeuvre par Talgorithme de type MAP unidirectionnel de la figure 2, 
30 appliquS k la detection de symboles dans un alphabet de A mots, ot 
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A = 2^; 

la figure 4 illustre les difKrentes 6tapes de calcul des variables "forward" 
mises en oeuvre par un algorithme de type MAP bidirectionnel selon 
rinvention ; 

5 - la figure 5 d6crit les diff^rentes 6tapes de calcul des variables "backward" 
mises en oeuvre par Talgorithme de la figure 3, ainsi que les diff6rentes 
Stapes de calcul des rapports de vraisemblance . 

Le principe g6n6ral de T invention repose sur la r^uction du nombre 
d'6tats s61ectionn6s lors du calcul du rapport de vraisemblance de Talgorithme 
10 MAP, de manifere k ^laborer un algorithme MAP simplifi6. Une telle induction est 
notanmient r6alis6e, dans un mode de realisation pr^fere de I'invention, par mise 
en oeuvre d'un algorithme de type "breadth-first", 

Prdalablement h la description d6taill6e des figures mentionn6es ci-dessus, 
on pr6sente le principe de la reduction d'6tats mise en oeuvre par Tinvention, ainsi 
15 que les diff6rentes equations Stabiles par les inventeurs. 

Dans un mode de realisation pr^fere de Tinvention, les inventeurs ont 
envisage d'avoir recours aux algorithmes de type "breadth-first", et notamment 
aux algorithmes de type M et de type T pour impl^menter des algorithmes MAP k 
nombre d'^tats r^duit 

20 Tout comme Talgorithme de Viterbi (d^crit notamment dans Farticle de G. 

D FORNEY « The Viterbi algorithm », Proceedings of the ffiEE, vol 61, No 3, 
march 1973), I'algorithme M est un algorithme de recherche en treillis de type 
« breadth-first » (en fitmgais "largeur d*abord"; un tel algorithme examine le cofit 
associe a une famille d'dtats k chaque instant) . L'algorithme de type M est 

25 notamment d6crit par J,B Anderson et S. Mohan dans « Sequential Coding 
algorithms : A survey and cost analysis » (en firangais, "Algorithmes de codage 
s6quentiels : enquSte et analyse de coflt"), IEEE Trans, on Conmiunications, vol. 
COM-32, Fevrier 1984. 
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Pour chaque profondeur du treillis de la figure 1, ralgorithme de Viterbi ne 
conserve que le meilleur chemin passant par chaque nceud, iUustr6 en trait 6pais 
sur la figure 1. 

L'algorithme de type M conserve quant h lui les M meilleurs chemins 
5 passant par les M meilleurs noeuds, oil M est un entier inf6rieur au nombre de 
nceuds du treillis de Talgorithme de Viterbi. Ainsi, lots de la mise en ceuvre de 
ralgorithme M, il est n^cessaire de classer les nceuds du treillis (symbolises par 
des points noirs sur le treillis de la figure 1) en fonction de leur m6trique cumul6e 
(on rappelle que la mStrique cumulfe des noeuds du treillis est 6quivalente k la 
10 variable "forward")* afin de s61ectionner les M meilleurs nceuds. 

L*algoridime de type T est 6galement un algorithme de type "breadth-first", 
mais diff&re de Talgorithme M par la m^thode mise en ceuvre pour s61ectionner 
les nceuds du treillis de la figure 1 propag6s par Talgorithme. En effet, 
ralgorithme T s^lectionne uniquement les nceuds du treillis associ6s k une 
15 m^trique cumul6e dont la distance k la meilleure m6trique cumul6e du treillis est 
infi^rieure k un seuil T. 

Selon rinvention, on prevoit avantageusement d^ntroduire de nouveaux 
algoritiunes T de type "breadth-first", k seuil variable. Pour ce faire, les inventeurs 
ont envisage deux approches distinctes : 
20 - dans une premifere variante de r6alisation de I'invention, on affecte un seuil, 
calculi a priori, k chaque transition (on rappelle qu\me transition d*un 6tat k 
un autre est symbolis6e par une branche du treillis de la figure 1). Ce seuil 
est calcul6 en tirant partie de la longueur finie de la sequence d' observation 
(sur le treillis de la figure 1, la sequence d'observation commence k I'instant 
25 k et se termine k Tinstant k+4). Ainsi, si Ton consid^re un treillis de 

longueur K+1, pour lequel BMAX repr^sente la valeur maximale des 
m^triques de branche, on choisit le semi T(k) tel que : ■ 
T(k)^(K+l-k)BMAX et 0<k<K ; 

dans une deuxi&me variante de realisation de Tinvention, on fixe la valeur du 
30 seuil de I'algorithme T, en mettant en ceuvre un algorithme adaptatif . En 
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effet, un avantage de ralgorithme de type T, par rapport h I'algorifhme M, 
est que, selon cet algorithme, il n'est pas n^cessaire de classer les noeuds du 
treillis en fonction de leur m^trique cumul^e, pour s^lectionner les meilleurs 
noeuds. L'algorithme de type T impose seulement de trouver le meilleur 
5 noeud du treillis, afin de nonnaliser les m6triques cumul6es de Tensemble 

des nceuds du treillis. Dans cette deuxifeme variante de realisation de 
rinvention, les inventeurs ont envisage d'utiliser un algorithme adaptatif 
pour fixer le seuil de Talgoritfame T, afin d^assendr sa complexity. 
Par exemple, en utilisant un algorithme de gradient, on pent mettre le seuil k 

10 jour par T^quation suivante : 

Tik + l)--nk)-fi(Mik)-Mc) (9) 
ot Mc repr6sente le nombre de noeuds propag6s vis6, Af(A:)le nomhre de noeuds 
propages k I'instant k, et oil A* est une petite constante positive que Ton appelle le 
gain d'apprentissage. 

15 On pent 6galement generaliser ralgorithme au cas oil depend de 

comme c'est le cas pour les algorithmes de type gradient k pas variable. 

On pr^sente desonnais plus en details la mise en oeuvre de tels algorithmes 
de type "breadth-first" dans le cadre de la reduction du nombre d'^tats des 
algorithmes de type MAP. 

20 Afin de pouvoir utiliser les algorithmes de type MAP avec un nombre 

d*6tats r^duit, Tinvention propose avantageusement d'affecter k chaque instant k 
une valeur par d6faut a(k) et b(Jc) respectivement aux variables "forward" et 

"backward", 

Solent m( et les ensembles des 6tats s61ectionn6s respectivement 

25 dans le sens "forward" et "backward" k Tinstant k. 

Si r^tat = i fait partie de Tensemble m( , c'est que cet 6tat a 6t6 

s61ectionn6 k Tinstant k. 

Les variables "forward" 6tant associ^es k des probability, elles doivent 

Verifier la condition de normalisation : 2«i* = 1- On peut d6montrer que si Ton 
30 veut verifier cette condition, alors il faut choisir a(k) tel que : 
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n apparalt de manifere 6vidente que les variables "forward" pour / g M/ne sont 
pas accessibles (les 6tats associds n'^tant pas s61ectionn6s, par definition de M( ), 

et Ton ne peut done que determiner des bomes de Tintervalle de variation de 
5 a(k) . En faisant Thypothfese que les Af 6tats 616ments de Af/ ont les M plus 

grandes probability a posteriori (ou APP), on borne a(k) par : 

0^a(Jk)^Mn(af) (11) 

Si Ton fixe a{k) k 0 quel que soit k , cela revient k consid&er que les 6tats 
non propag6s ont une APP nulle et ne participent pas au calcul de la 
10 vraisemblance donn6 par T^quation (8). En d'autres termes, cela revient k attribuer 
une confiance absolue aux choix des M 6tats propag6s. 

A contrario, si Ton fixe a(k) k sa valeur maximale, alors on attribue une 

confiance minimale au choix des M 6tats propag^ . 

On peut s^pliquer le m8me raisonnement aux variables "backward", ce qui 
15 conduit k determiner b(k) tel que : 

0<b{k)<Min{P^) (12) 

Les equations (11) et (12) ci-dessus peuvent 6tre affinees dans le cadre de 
Tutilisation d'un algorithme de type M d'une part, et d'un algorithme k seuil de 
type T d'autre part. 

20 En effet, lorsque Ton met en oeuvre un algorithme de type M, on dispose k 

chaque instant jt, du classement des etats seiectionn6s dans le sens direct (ou 

"forward") et indirect (ou "backward") respectivement, en fonction des valeurs de 
a * et fi* . n est done aise de determiner les valeurs de Mm(af ) et Mm(i3f ), et 

done d'affecter k a(k) et b(k) les valeurs suivantes : 
25 aik) = Mm{af)-Cj. (13) 

K*) = M/Ai(i8*)"C^ (14) 

avec Cy.etCj deux constantes positives. Ces deux constantes Cf et permettent de 
regler la confiance que Ton accorde au choix des etats propages. 
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Dans le cas d'un algorithme de type T, c'est-k-dire lorsque les 6tats 
propag6s par Talgorithme sont s61ectionn€s par un seuil (fixe ou variant au cours 
du temps conime propos6 pr6c6deimnent), la dfitermination des minima des 
Equations (11) et (12) necessite une charge de calcul suppl6mentaire importante. 
5 Par consequent, les inventeurs proposent, pour cette classe d' algorithme, de fixer 
a(^)et b(k) par rinterm6diaire du seuil : 

a{k)^T^ik)-Cj (15) 
bik)^T\k)-c^ (16) 
oil (it) reprfeente le seuil dans le sens "forward" k r instant k, et oii T^(k) 

10 repr6sente le seuil dans le sens "backward" & Tinstant k (tels que d^finis par 
exemple par F^uation (9) ci-dessus), et Cy 6tant deux constantes positives. 

Dans la pratique, pour des raisons de simplicity d'impl6mentation, les 
algorithmes de type Log-MAP ou Max-Log-MAP sont souvent pr6f6r& aux 
algorithmes MAP classiques. En effet, le calcul du rapport de vraisemblance tel 

15 que d^fini par T^quation (3) ci-dessus necessite, pour les algorithmes de type 
MAP classiques, d'effectuer des divisions qui sont couteuses et peuvent poser des 
probldmes de stability. Pour cette raison, on pr6f6re g6n6ralement utiliser les 
algorithmes travaillant dans le domaine logarithmique, qui propagent Toppos^ du 
logarithme des variables "forward" et "backward", et aboutissent au calcul du 

20 logarithme du rapport de vraisemblance, encore appel6 Log-vraisemblance (ou 



LLR pour "Log-Iikelihood-Ratio") d6fini par : 



(3bis) 



oil PA/5? = P[X, = iVFo = yo^-'^^K = ,0 ^ ^ ^ . 

Si Ton consid^re un codeur convolutif g6n6rant des symboles X^, transmis 
25 par un canal bruit6 vers un r6cepteur qui regoit des symboles et impl6mente un 
algorithme de type Lx)g-MAP pour decoder les bits regus, le signe de LU^ donne 

la valeur du bit d6cod6, et sa valeur absolue indique la quality de la decision. 
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L'application du pxincipe de la reduction d'6tats, pr6sent6 ci-dessus pour 
les algorithmes MAP classiques, aux algoiithmes de type Log-MAP ou Max-Log- 
MAP est pr^ent^e bii^vement ci-apr6s. 

Solent af =-Log(ocf) et les variables logarithmiques 

5 "forward" et "backward". Selon I'invention, pour les algorithmes de types Log- 
MAP et Max-Log-MAP k nombre d'^tats r^duit, les variables logarithmiques 
manquantes sont remplacdes dans le sens "forward" et "backward" lespectivement 
par aik) et b(k) d^finis par : 

aik)^Max(a'^) (17) 

tsuf 

10 Hk)>Max(p',) (18) 

leAff 

Ainsi, selon rinvention, lorsque la reduction du nombre d*6tats de 

Talgorithme de type Log-MAP ou Max-Log-MAP est r6alis6e par mise en oeuvre 

d'un algorithme "breadth-first" utilisant un algorithme de tri (c'est-Ji-dire un 

algoritfame de type M), on remplace les donnas manquantes respectivement dans 

15 le sens "forward" et "backward" par : 

a(k)^MaxicCi)^Cf (19) 

ieMf 

b(k) = Max(p'!)-^ct (20) 

avec Cf et Cb deux constante positives. 

En ce qui conceme les algorithmes "breadth-first" utilisant un seuil pour 
20 s61ectionner les 6tats propag^s dans le sens "forward" et "backward" (c'est-k-dire 
les algorithmes de type T), les donn^es manquantes sont remplac6es 
respectivement par : 

a(fc) = r^(fe) + c/ (21) 
Kife) = T*(fe) + C6 . (22) 

25 avec T^(k) et (it) respectivement les seuils "forward" et "backward" (variant 
6ventuellement en fonction du temps k , conome propos6 pr6c6demment). 

^invention permet 6galement de combiner au principe de la reduction 
d*6tats pr^sent^ pr6c6denunent» des aspects relatifs k la simplification des 
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algoriihmes classiques de type MAP (c'est-i-dire MAP, Log-Map et Max-Log- 
MAP), 

Comme mentioim6 pr6c6demment, les algorithmes MAP, Log-MAP et 
Max-Log-MAP n6cessitent la propagation des variables "forward" et "backward" 
5 respectivement dans le sens direct et indirect, par rapport au treillis associ6. 

Dans une variante de realisation de Tinvention, on propose de simplifier 
ces algorithmes en fixant les variables "forward" ou "backward" h une valeur 
pr6d6ternun6e, de pr6f6rence constante, de fa5on k elaborer des algorithmes dits 
unidirectionnels. Ainsi, en fixant les variables "forward", a* ou af k une valeur 

10 constante quelconque quels que soient k et on obtient des algorithmes MAP, 
Log-MAP et Max-Log-MAP dits "indirects". De meme, en fixant la variable 
"backward", jS* ou k une valeur constante quelconque quels que soient k et 

/ , on obtient des algorithmes dits "directs". 

On peut alors appliquer le principe de la reduction d*6tats k ces algorithmes 

15 directs et/ou indirects, en mettant en ceuvre les algorithmes de type "breadth-first" 
d^crits pr6c6demment, de fa^on k obtenir des algorithmes directs et/ou indirects k 
nombre d'6tats r6duit. Ainsi, dans le cas d'un algorithme MAP indirect par 
exemple, pour lequel toutes les variables "forward" af ont 6t6 fix6es k une valeur 
constante, on remplace les variables '"backward" manquantes jSf (correspondant k 

20 des etats non s61ectionn6s) par une valeur b(k) pr6detenmn6e, ainsi que d6crit 
pr6c6demment dans le document 

De telles m6thodes de gestion des donn^es manquantes permettent de 
rendre efficaces les algorithmes directs et indirects utilisant un nombre d'€tats 
r6duit. En effet, lorsque le nombre d'6tats propag6s est faible, il est probable que 

25 Tun au moins des ensembles sur lequel on calcule le rapport de vraisemblance 
d^fini par T^quation (3) soit vide. 

La m^thode propos6e, selon laquelle on affecte des valeurs par d^faut aux 
variables intermddiaires n^essaires au calcul du rapport de vraisemblance mais 
correspondant k des 6tats non propag^ du treillis, permet alors de determiner le 

30 rapport de vraisemblance d&s qu'il y a au moins deux dtats propag6s dans le 
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treillis, et de lui affecter une valeur souple (encore appel^ valeur "soft"). En effet, 
en affectant la valeur 0 k la probability d'un 6tat non propag6, on obtient, pour les 
algorithmes de type MAP, des valeurs de rapport de vraisemblance en ^^jj^^^ 
coDStenta^ pour les algorithmes de type Log-MAP, on obtient des valeurs 

5 de rapport de vraisemblance en — oo ou en +«> . 

Un inconvenient de ces algorithmes unidirectionnels, que les inventeurs 
ont identifi6, est que les decisions prises par de tels algorithmes, c'est-k-dire le 
signe des rapports de vraisemblance d6finis par Tfiquation (3) (ou (3bis)), different 
g6n6ralement des decisions prises par les algorithmes mettant en oeuvre un "trace- 

10 back". (On rappelle que les algorithmes mettant en oeuvre un "trace-back" 
permettent une prise de d^ision correspondant au chemin le plus vraisemblable). 

En effet, les 6tats associ^s aux decisions prises par les diff^rents algorithmes 
unidirectionnels peuvent ne pas §tre autoris6s par la matrice de transition 11 . Dans 
une variante de realisation avantageuse de Tinvention, on a done envisage 

15 d'utiliser conjointement un algorithme unidirectioimel k nombre d*etats r^duit 
avec un "trace back" comparatif ou comparatif pond^re* 

Selon rinvention, le "trace back" comparatif consiste k comparer les 
decisions prises par un algorithme unidirectionnel aux decisions correspondantes 
prises par un algorithme de type Viterbi, ainsi qu'illustre en figure 2. On considdxe 

20 y^xm ensemble d'etats associes k la decision D[ foumie par Talgoritfune de type 
Viterbi k I'instant L Soil A^, la vraisemblance associee k ^ k Tmstant ^, 
calcuiee par Talgorithme unidirectionnel et la decision correspondante prise 

par Valgorithme unidirectionnel. 

Au cours d'une premiere 6tape referencee 21, on calcule le rapport de 
25 vraisemblance A^, de fagon k prendre une decision A^ pour Talgorithme 
unidirectionnel. On analyse (22) la decision Df correspondante pour Talgorithme 
de type Viterbi et on compare (23) ces deux decisions A^ et D[ • 

Si les deux decisions dif^rent, alors la decision prise par Talgorithme de 
type Viterbi remplace (25) celle de Talgorithme unidirectionneL La valeur absolue 
30 des rapports de vraisemblance associes k ces decisions reste inchangee ou forcee k 
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une valeur V . Vest soit fix6e a priori k une valeur petite soit calcul6e pour 
chaque sequence, Dans ce dernier cas, on pent par exemple fixer la valeur 
absolue du rapport de vraisemblance moyen de la s^uence. 

Plutdt que de choisir brutalement entre la d&:ision foumie par Talgoritiune 
5 de type Viterbi Df et celle produite par I'algorithme unidirectionnel A^, on pent 

6galement faire une pond^ration de ces deux decisions : c'est ce que Ton appelle le 
principe du "trace back" comparatif pond6r6, 6galement illustr6 en figure 2. 

Lors de T^tape de "trace back" comparatif pond6r6, la vraisemblance est 
remplac6e (24) par A^ d6finie par : 
10 Ay=A^ + A''xV (23) 

avec V 6gale a la valeur absolue du rapport de vraisemblance moyen calculee sur 
la sequence, ou Sx6o a priori. 

On pr^sente dSsormais, en relation avec la figure 3, un mode de r&disation 
des diff6rentes Stapes de calcul mises en oeuvre par un algorithme de type MAP 
15 logarithnaique unidirectionnel selon Tinvention, appliqu6 a la detection de 
symboles dans un alphabet de A mots, oft A = 2^. 

Le vecteur CM(i) 31 contient les M(i) m^triques cumul6es correspondant 
aux M(i) etats retenus a F instant i stock6s dans le vecteur State(i) 32 de taille M(i) 
foisQbits. 

20 Le bloc 33 « calcul BM et ACS » calcule les A.M(i) mdtriques cumul6es 

assocites aux A.M(i) 6tats suivant les M(i) 6tats du vecteur State(i) 32. n efifectue 
ensuite la selection des meilleurs 6tats convergeant vers le mSme noeud. Cette 
selection consiste soit h s^lectioimer les 6tats qui ont la plus petite m^trique 
cumulee (dans le cas d'un algorithme de type Max-log-map unidirectionnel) soit k 

25 effectuer r6cursivement une operation de comparaison correction en utilisant la 
formule du logarithme Jacobien (dans le cas d'un algorithme de type Log-Map 
unidirectionnel). 

Ainsi, dans le cas d'un algorithme unidirectionnel de type Max-Log-MAP, 
on procMe en trois Stapes successives : 
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a) Premifere 6tape : calcul de la m6trique de branche 
BM(ij>^ln(fej(y*)), conform^ment k T^quation (5) d6crite pr6c6demment On 
notera que la metrique de branche s*exprime simplement comme une distance 
euclidienne. 

5 b) Deuxieme 6tape : accumulation. 

On met k jour les variables (respectivement P^) de manidre simplifi^e. 
On considere, h, Tinstant k-1, deux noeuds Pq et Pi associ6s respectivement aux 
variables "forward" aj^ et a^. A Tinstant k suivant, on considfere le nceud final 

associ6 k la m6trique cumul6e CMj(k). La m6trique de branche du nceud Pq 

10 (respectivement P,) au nceud final est BM(Po(j)d) (respectivement BM(Pi(j) j). On 
a alors : ai^^Jisa{^^ 

c) Troisifeme dtape : selection. 

On choisit le meilleur chemin qui converge vers chaque nceud : c'est cette 
6tape qui va permettre le "trace-back". 
15 Pour un algorithme de type Log-MAP, les premiere et troisi^me Stapes 

sont identiques, mais, au cours de la deuxieme 6tape de mise k jour des variables 
"forward", on a : 

oil f(z)=ln(l+exp(-z)). 

20 Les m6triques des meilleurs noeuds s^ectionn^s par le bloc 33 « calcul BM 

et ACS » sont sauvegard6es dans le vecteur TmpCM(i+l) 34. Le nombre d'^tats 
s61ectionn6s k cette 6tape est compris entre M(i) et A-M(i). Les 6tats 
correspondants s61ectionn6s k cette etape sont stock^s dans le vecteur 
TmpState(i+l) 35. 

25 Le bloc 36 « selection M ou T » pemiet de s61ectionner les 6tats propag^s 

par Talgorithme MAP. 

Si la reduction d'6tats met en ceuvre un algorithme de type T, le bloc 36 
« selection M ou T » s61ectionne les £tats dont la distance de la metrique cumul6e 
par rapport k la meilleure m6trique cumul6e est infi^rieure k un seuil T. 
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Si la reduction d'6tats met en oeuvre un algoiithme de type M, le bloc 36 
« selection M ou T » s61ectionne les M meilleurs 6tats, 

Dans le cas de ralgorithme de type T, le bloc 36 « selection M ou T » 
calcule 6galement le nouveau seuil T(i+1) en fonction du nombre d'6tats 
5 s61ectionn6s M(i) en appUquant le seuil T(i), comme d6crit pr6c6deimnent dans 
r^quation (9). 

A Tissue de T^tape de selection effectu^e par le bloc r6f6renc6 36, on 
stocke dans le vecteur CM(i+l) 37, les M(i+1) metriques cumul^es correspondant 
aux M(i+1) 6tats retenus k I'instant i+1 stock6s dans le vecteur State(i+1) 38 de 

10 taille M(i+1) fois Q bits. 

Le vecteur Survivor(i+l) 39, de taille N mots de Q bits contient Tindice 
des M(i) branches qui convergent sur les M(i) 6tats retenus, stock6s dans le 
vecteur State(i) 32. L'6tat sur lequel converge la branche gagnante du treillis est 
cod6 par Tadresse de stockage de Tindice de la branche, sur Q bits. 

15 C'est avec Tensemble des vecteurs Survivor(i+l) 39, pour i variant de 0 k 

K-1, que Von effectue I'opdration de "trace back", afin de determiner les decisions 
k comparer ou a ponderer avec celles determindes par le bloc 40 de calcul du ULR 
(pour "Log-Likelihood Ratio", ou rapport de vraisemblance logarithmique), ainsi 
que d6crit prec^demment en relation avec la figure 2. 

20 On pr6sente d6sormais, en relation avec la figure 4, un mode de realisation 

du calcul des variables "forward" pour un algorithme logarithmique MAP 
bidirectionnel selon Tinvention. 

On rappelle que les algorifhmes bidirectionnels de type Log-Map ou Max- 
Log-Map k nombre d'6tats r^duit commencent, dans une premidre 6tape, par 

25 calculer r^cursivement les variables "forward", comme pour les algorithmes 
unidirectionnels. 

La figure 4 illustre le calcul de la liste des 6tats conserves et des metriques 
cumulees dans le sens direct ou "forward". 

A la differences des algorithmes unidirectionnels pr6sent6s en figure 3, les 
30 m6triques cumulees du vecteur CMF(i) 42, les etats s^lectionnes correspondants 
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du vecteur StateF(i) 43, et 6ventuellement le nombre d*6tats s61ec1ioim6s k chaque 
instant, sont stock^s, afin d*Stre utilises pour calculer les rapports de 
vraisemblance LLR lors de r6tape "backward". 

Comme iUustre par la figure 3, le vecteur CMF(i) 42 contient les MF(i) 
5 metriques cumul6es correspondant aux MF(i) 6tats retenus dans le sens "forward" 
k rinstant i stock6s dans le vecteur StateF(i) 43 de taille MF(i) fois Q bits. 

Le bloc 33 « calcul BM et ACS » calcule les A.MF(i) metriques cumul6es 
associ^es aux A.MF(i) 6tats suivant les MF(i) 6tats du vecteur StateF(i) 43. n 
effectue ensuite la selection des meilleurs 6tats convergeant vers le m8me nceud, 
10 selon un proc6d6 similaire k celui d6crit en relation avec la figure 3. 

Les m6triques des meilleurs noeuds s61ectionn6s par le bloc 33 « calcul BM 
et ACS » sont sauvegard^s dans le vecteur TmpCMF(i+l) 44. Le nombre d*6tats 
s61ectionnes k cette 6tape est compris entre MF(i) et A.MF(i). Les 6tats 
correspondants s61ectionn6s k cette 6tape sont stock^s dans le vecteur 
15 TmpStateF(i+l) 45. 

Le bloc 36 « selection M ou T » permet de s61ectionner les etats propag^s 
par Talgorithme MAP. 

Si la reduction d'6tats met en oeuvre un algorithme de type T, le bloc 36 
« selection M ou T » s61ectionne les 6tats dont la distance de la m^trique cumul6e 
20 par rapport k la meiUeure m^tcique cumulde est inf6rieure k un seuil T. 

Si la reduction d'6tats met en oeuvre un algorithme de type M, le bloc 36 
« s^ection M ou T » s61ectionne les M meilleurs etats. 

Dans le cas de ralgorithme de type T, le bloc 36 « selection M ou T » 
calcule 6galement le nouveau seuil dans le sens "forward" TF(i+l) en fonction du 
25 nombre d'6tats s61ectionn6s MF(i) en appliquant le seuil TF(i), comme d6crit 
pr6cedemment dans I'equation (9). 

A Tissue de r6tape de selection effectu6e par le bloc r6f6renc6 36, on 
stocke dans le vecteur CMF(i+l) 46, les MF(i+l) m6triques cumul6es 
correspondant aux MF(i+l) 6tats retenus k Tinstant i+1 stock^s dans le vecteur 
30 StateF(i+l) 47 de taiUe MF(i+l) fois Q bits. 
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Le calcul des m6triques cumul^es **backwaid" du vecteur CMB(i) SO et des 
6tats s61ectioim6s stock6s dans le vecteur StateB(i) 51 s'effectue de mani^re 
sym6trique au calcul "forward", ainsi qu'illustr6 par la figure 5. La figure 5 
presente 6galement le calcul des LLR h, partir des variables "forward" (mdtriques 
5 cumul6es, liste d'6tats et 6ventuellement seuils) m6morisfes dans le sens direct, et 
des nouvelles variables calculees dans les sens indirect 

Ainsi, le vecteur CMB(i) 50 contient les MB(i) m6triques cumulees 
correspondant aux MB(i) 6tats retenus dans le sens "backward" k Tinstant i 
stock6s dans le vecteur StateB(i) 51 de taille MB(i) fois Q bits. 
10 Le bloc 33 « calcul BM et ACS » calcule les A.MB(i) m6triques cumuli 

assocites aux A.MB(i) 6tats suivant les MB(i) 6tats du vecteur StateB(i) 51. n 
effectue ensuite la selection des meillemrs 6tats convergeant vers le mSme nceud, 
selon un proc6d^ similaire k celui ddcrit en relation avec les figures 3 et 4. 

Les mStriques des meilleurs noeuds sSlectionn^s par le bloc 33 « calcul BM 
15 et ACS » sont sauvegard6es dans le vecteur TmpCMB(i+l) 53. Le nombre d'^tats 
s61ectionn6s h cette 6tape est compris entre MB(i) et A.MB(i). Les 6tats 
correspondants sdlectionn6s k cette 6tape sont stock^s dans le vecteur 
TmpStateB(i+l) 52. 

Le bloc 36 « selection M ou T » permet de s61ectionner les 6tats propag^s 
20 par Talgorithme MAP. 

Si la reduction d'6tats met en oeuvre un algorithme de type T, le bloc 36 
« selection M ou T » s61ectionne les 6tats dont la distance de la m6trique cumul6e 
"backward" par rapport k la meilleure m^trique cumuli est infdrieure k un seuil 
T. 

25 Si la reduction d'6tats met en oeuvre im algorithme de type M, le bloc 36 

« selection M ou T » s61ectionne les M meilleurs 6tats. 

Dans le cas de Talgorithme de type T, le bloc 36 « selection M ou T » 
calcule 6galement le nouveau seuil dans le sens "backward" TB(i+l) en fonction 
du nombre d'^tats s^lectionn^s MB(i) en q>pliquant le seuil TB(i), comme d6crit 

30 pr&:&lemment dans T^quation (9). 
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A Tissue de T^tape de selection effectu^e par le bloc T6f6renc€ 36, on 
stocke dans le vecteur CMB(i+l) 54, les MB(i+l) m^triques cumul^es 
correspondant aux MB(i+l) 6tats retenus k Tinstant i+1 stock6s dans le vecteur 
StateB(i+l) 55 de taille MB(i+l) fois Q bits. 
5 Le calcul du LLR 40 utUise les m^triques cumulfes du vecteur CMF(i) 42 

et du vecteur CMB(i) 50, ainsi que les 6tats s61ectionn6s stock^s dans les vecteurs 
StateF(i) 43 et StateB(i) 51. Dans le sens "forward" les m^triques cumul^es du 
vecteur CMF(i) 42, les 6tats du vecteur StateF(i) 43, et 6ventuellement les seuils 
TF(i) et le nombre d'6tats s61ectionn6s MF(i), sont m6moris6s pour toute la 

10 sequence, c*est-^-diie pour i variant de 1 ^ K. 

Le bloc « Calcul du LLR » 40 pent utiliser les seuils TF(i) et TB(i), dans le 
cas d'une reduction d'6tats par algorithme de type T, ou les moins bonnes 
m^triques cumul6es CMB(M) et CMF(M), dans le cas d*une reduction d'^tats par 
algorithme de type M, pour g6rer les donn6es manquantes, comme expose 

15 pr6c6deinment. Tout comme pour les algorithmes unidirectionnels, le calcul du 
LLR utilise T approximation Max-Log-Map ou la formule du logarithme Jacobien 
pour les algorithmes de t^pe Log-Map. 
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REVENDICATIONS 

1. Proc6d6 de traitement d'un signal mettant en oeuvre un algorifhme de type 
MAP ("Maximum a posteriori") permettant de determiner un rapport de 
vraisemblance A^^ d'un ensemble d'6tats X d'un treillis h un instant k, 

5 k chacun desdits 6tats 6tant associ6e au moins une variable interm6diaire, 
appartenant au groupe comprenant une variable appel^ "forward" et une variable 
appel6e "backward", propag^es par ledit algorithme MAP et calcul6es 
r^cursivement respectivement dans un sens direct et dans un sens indirect audit 
instant k par rapport audit treillis, 

10 caract6ris6 en ce qu'U comprend une 6tape de reduction du nombre d'6tats 
s61ectionn6s par ledit algorithme de type MAP en vue d'lm calcul dudit rapport de 
vraisemblance, 

et en ce que, pour au moins certains Stats non-s€lectionn6s, on affecte h ladite 
variable "forward" et/ou "backward" correspondante au moins une valeur 
15 d6termin6e, de fagon k calculer un rapport de vraisemblance approchS. 

2. Proced6 selon la revendication 1, caracteris6 en ce que, pour un instant k 

donn6, ladite au moins une valeur d6tennin6e a(k) affectSe h ladite variable 
"forward" est telle que 0<a(k)^Min(ah, et/ou ladite au moins une valeur 

d6termin6e b(k) affectee k ladite variable "backward" est telle que 
20 0^fc(Jk)<Mm(A*). ovL M/ et Mk'' reprfisentent un ensemble desdits 6tats 

sglectionnSs respectivement dans ledit sens direct et dans ledit sens indirect audit 
instant k, et oil al" et repr6sentent respectivement lesdites variables "forward" 
et "backward" audit instant k. 

3. Proc6d6 selon la revendication 2, caract6ris6 en ce que, pour un instant k 
25 donn6, ladite valeur d^terminee a(k) et/ou b(k) est unique et est affectee k au 

moins une variable "forward" Oi*" et/ou "backward" Pi^ 

4. Proc6d6 selon Tune quelconque des revendications 1 k 3, caract6ris6 en ce 
qu'on affecte une valeur constante auxdites variables "forward", respectivement 
"backward", de fagon que ledit algorithme de type MAP soit un algorithme 

30 unidirectionnel de type direct, respectivement indirect 
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5. Proc6d6 selon Tune quelconque des revendications 1 k 4, caract6ris6 en ce 
que ladite 6tape de reduction du nombre d'^tats met en oeuvre un algoiithme de 
recherche en treillis de type "breadth-first", 

6. Proc6d6 selon la revendication 5, caract6ris6 en ce que ledit algorithme de 
5 type "breadth-first" est un algorithme de type M. 

7. Proc6d6 selon la revendication 5, caract6ris6 en ce que ledit algorithme de 
type "breadth-first" est un algorithme de type T, utilisant au moins un seuil. 

8. Proc6d6 selon la revendication 7, caract6ris6 en ce que ledit au moins un 
seuil est variable en fonction dudit instant k. 

10 9. Proc6d6 selon la revendication 8, caract6ris6 en ce qu*on affecte audit seuil 
variable une valeur prdd6tennin6e pour chaque instant k. 

10. Proc6d6 selon la revendication 8, caract6ris6 en ce que, pour chaque 
instant k, la valeur dudit seuil variable est d6tenmn6e par mise en ceuvre d'un 
algorithme adaptatif . 

15 11. Proc6d6 selon la revendication 10, caract6ris6 en ce que ledit algorithme 
adaptatif est du tj^e algorithme de gradient. 

12. Proc^de selon Tune quelconque des revendications 10 et 11, caract6ris6 en 
ce que, ledit treillis comprenant une plurality de noeuds associes chacun k un 
desdits 6tats et k un instant k donn6, la valeur dudit seuil variable T £t un instant 

20 (k+1) est d6temiin6e par I'^quation suivante : 

nk + 1) = nk) - iimk) - m,) 

oii T(k) repr6sente la valeur dudit seuil variable audit instant k. 
Me est le nombre de noeuds dudit treillis propag6s vis6, 
MCk) est le nombre de noeuds dudit treillis propag6s audit instant k, 
25 et II est une constante positive representative dhm gain d^apprentissage. 

13. Proc6d6 selon Tune quelconque des revendications 11 et 12, caract^ris^ en 
ce que ledit algorithme adaptatif est du type algorithme de gradient h pas variable. 

14. Proc6d6 selon Tune quelconque des revendications 12 et 13, caract6ris6 en 
ce que ledit gain d'apprentissage \i est fonction dudit instant k. 
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15. Proc6d6 selon Tune quelconque des revendications 2 k 14, caract6ris6 en ce 

que, ledit algorithme de type "breadth-first" 6tant un algorithme de type M, 

lesdites valeurs d6tenmn6es a(k) et/ou b(k) affect^es respectivement auxdites 

variables "forward" et/ou "backward", k un instant k donn6 sont donn6es par les 

5 Equations suivantes : 
a(k) = Min{ah - c. 

bik) = Mi]^iP^)-c, 

oh Cf et Cb sont deux constantes positives. 

16. Proc6d6 selon Tune quelconque des revendications 2 k 14, caract6ris6 en ce 
10 que, ledit algorithme de type "breadth-first" dtant un algorithme de type T, lesdites 

valeurs d6temain6es a(k) et/ou b(k) afifect^es respectivement auxdites variables 

"forward" et/ou "backward", k un instant k donn6, sont donn^ par les Equations 

suivantes : 

a(k) = T^ik)-Cj. 

15 b{k)^T^(k)-c^ 

ot Cf et Cb sont deux constantes positives, et ot T^(k) et T^Ck) designent la valeur 

dudit seuil variable audit instant k respectivement dans ledit sens direct et dans 

ledit sens indirect. 

17. Proc6d6 selon Tune quelconque des revendications 1 k 16, caract6ris6 en ce 
20 que ledit algorithme de type MAP appartient au groupe comprenant : 

les algorithmes de type MAP ; 
les algorithmes de type Log-MAP ; 
les algorithmes de type Max-Log-MAP. 

18. Proc6d6 selon Tune quelconque des revendications 4^17, caract6ris6 en ce 
25 que, ledit algorithme de type MAP 6tant un algorithme unidirectionnel, ledit 

proc6d6 met en oeuvre une ^tape de comparaison de decisions prises par ledit 
algorithme unidirectionnel avec des decisions correspondantes prises par un 
algorithme de type Viterbi, appel^es d^sions de Viterbi. 

19. Proc6d6 selon la revendication 18, caract6ris6 en ce qu'en cas de 
30 comparaison negative pour au moins une desdites decisions prises par ledit 
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algorithme unidirectioimel, U met en ceuvre une 6tape de substitution de ladite 
decision de Viterbi correspondante k ladite decision prise par ledit algorithme 
unidirectionnel, appel6e decision substitute. 

20. Proctde selon la revendication 19, caract6ris6 en ce qu'on affecte k la 
5 valeur absolue dudit rapport de vraisemblance associt k ladite d6cision substitute 

une valeur dtterminee V. 

21. Procddt selon la revendication 20, caracttrist en ce que ladite valeur 
determinee V est tgale k la valeur absolue du rapport de vraisemblance moyen de 
la sequence. 

10 22. Proctde selon la revendication 18, caracttrist en ce qu'en cas de 
comparaison negative pour au moins une desdites decisions prises par ledit 
algorithme unidirectionnel, il met en ceuvre une 6tape de pondtration dudit 
rapport de vraisemblance associt k ladite decision considtrte, tenant compte de 
ladite decision de Viterbi. 

15 23. Procede selon la revendication 22, caracttrist en ce que, Y ttant un 
ensemble d'ttats associts k une decision foumie par ledit algorithme de type 
Viterbi a un instant i, et representant le rapport de vraisemblaace associt kYk 
rinstant i, calcul6 par ledit algorithme unidirectionnel, au cours de ladite 6tape de 
pondtration, on remplace A^^ par dtfini par = + x V, 

20 oh V est une valeur dtterminte. 

24. Application du proc6d6 selon Tune quelconque des revendications Ik 23 k 
Tun au moins des domaines appartenant au groupe comprenant : 

la detection de symboles ; 
le codage/dtcodage de signaux ; 
25 - le turbo-dtcodage ; 

la turbo-dttection ; 

le codage de source par quantification en treillis. 

25. Recepteur de signaux de communication comprenant des moyens de mise 
en oeuvre d'un algorithme de type MAP ("Maximum a posteriori") pennettant de 
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d^temiiner un rapport de vraisemblance d'un ensemble d'^tats X d'un treillis k 
un instant k, 

k chacun desdits €tats dtant associ6e au moins une variable interm6diaire, 
appartenant au groupe comprenant une variable appel^e "forward" et une variable 
5 appel6e "backward", propag^es par ledit algorithme MAP et calcul6es 
r^cursivement respectivement dans un sens direct et dans un sens indirect audit 
instant k par rapport audit treillis, 

caract6ris6 en ce qu'il comprend des moyens de reduction du nombre d'6tats 
s61ectionn€s par ledit algorithme de type MAP en vue d'un calcul dudit rapport de 
10 vraisemblance, 

et en ce que, pour au moins certains 6tBts non-s61ectionn6s, on affecte k ladite 
variable "forward" et/ou "backward" correspondante au moins une valeur 
d^termin^, de f agon k calculer un rapport de vraisemblance approch6. 
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